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Context. The study of stellar multiple systems provides us with important information about the stellar formation processes and 
can help us to estimate the multiplicity fraction in the Galaxy. 65 UMa belongs to a rather small group of stellar systems of higher 
multiplicity, whose inner and outer orbits are well-known. This allows us to study the long-term stability and evolution of the orbits 
in these systems. 

Aims. We obtained new photometric and spectroscopic data that when combined with interferometric data enables us to analyze the 
system 65 UMa and determine its basic physical properties. 

Methods. We perform a combined analysis of the light and radial velocity curves, as well as the period variation by studying the times 
of the minima and the interferometric orbit. A disentangling technique is used to perform the spectra decomposition. This combined 
approach allows us to study the long-term period changes in the system for the first time, identifying the period variation due to the 
Ph motion on the visual orbit, in addition to some short-term modulation. 

Results. We find that the system contains one more component, hence we tread it as a sextuple hierarchical system. The most inner 
pair of components consists of an eclipsing binary orbiting around a barycenter on a circular orbit, both components being almost 
identical of spectral type about A7. This pair orbits on an eccentric orbit around a barycenter, and the third component orbits with 
a period of about 640 days. This motion is reflected in the period variation in the minima times of the eclipsing pair, as well as in 
the radial velocities of the primary, secondary, and tertiary components. Moreover, this system orbits around a barycenter with the 
distant component resolved interferometrically, whose period is of about 118 years. Two more distant components (4" and 63") are 
| also probably gravitationally bound to the system. The nodal period of the eclipsing-pair orbit is on the order of only a few centuries, 

which makes this system even more interesting for a future prospective detection of changing the depths of minima. 
Conclusions. We identify a unique solution of the system 65 UMa, decomposing the individual components and even shifting the 
system to higher multiplicity. The study of this kind of multiple can help us to understand the origin of stellar systems. Besides 
65 UMa, only another 1 1 sextuple systems have been studied. 

Key words, binaries: eclipsing - stars: fundamental parameters - stars: early-type - stars: individual: 65 UMa, DN UMa, HD 103795, 
I/-). SAO 43913 
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[~^! ■ 1. Introduction The study of systems of higher multiplicity is still rela- 
(-—•) tively undeveloped yet, and can provide insigh t into their for- 
; As members of more complex multiple systems, the eclipsing mation MoreoV er, iGoodwin & KroTu^l (£005) found that the 
binaries can provide us important information about their phys- majority of the early-type stars are found in multiple systems. 
.. . ical properties, as derived from different methods. This is the star-forming theories are still based on many ad hoc assump- 
m j> ; case for 65 UMa, a system whose the close components form an tions and the physica i characteristics of the multiple systems 
k> . eclipsing binary, and additional components foundjo be grav- can provide strong constraints on some of them . These can be 
"£j ■ nationally bounded to this pair ( |Pourbaix et alj. |2004j). Thanks e g the mass ratios of the inner and ou t er pairs, the ratio of p eri- 
i to the combined analysis, we have been able to derive the radii, ods , a nd in clinations, see for instance (Go odwin et al.ll2QQ7l and 
masses, and evolutionary statuses of the close components, in |Tokovininll2008h . In addition, the multiplicity fraction is one of 
addition to some properties of the distant ones. These systems the most cmcial param e ter s in theoretical models and nowadays 
are still very rare and mostly he relatively close to the solar sys- we kn ow of only 20 quintuple s, 1 1 sextuples, and 2 septuple sys- 
tem. Only 39 such systems are known where a close eclipsing temfj (Eggl eton & Tokovininll2008l) . 
binary is a member of a wide visual binary and we know both or- 
bits, their mutual inclinations, ratio of periods, etc. For instance, 

the ratio of periods can tell us something about the long-term 2 The System 65 UMa 

stability of the system. These unique systems are the most suit- 
able for studies of dynamical effe cts, the short and long-term The multiple system 65 UMa (= WDS Jl 1551+4629) consists 
evolution of the orbits, etc. (see e.g. Soderhie lmll975l) . of four visible components. The angular distance between the 

primary and component D (= HR 4561) is about 63", while 

* Tables Al to A4 are only available in electronic form at the the C component is at a distance of about 4". The primary 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 65 UMa AB = DN UMa, the brightest member of the sys- 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ tem, was also resolved to be a binary via classical micromet- 
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ric measurements by lAitkenI (Il908l) . Since then, many precise 
interferometric observations of this close pair were carried out. 
More over, the primary componen t was di scovered to be a vari- 
able (IGimenez & OuesadaLTl979l) . Later, iGimenez & Ouesa da 
( 1982j) found the star to be an eclipsing binary with the orbital 
period of about 1.7304 d. Its variability was not discovered ear- 
lier owing to its rather shallow eclipses of only about 0.09 mag. 
This is caused by the presence of other components in the aper- 
ture and therefore a large fraction of the third light, which re- 
duces the photometric amplitude of t he eclipses. The light curve 
(hereafter LC) analysis presented by Garci a & Gimenej {l986) 
was slightly indicative of an eccentric orb it. The radia l veloc- 
ity curve (hereafter RV) was analyzed bv lPopperl (Il986l) . All of 
these studies found that the close binary consists of two very 
similar stars. 

The position angle between the A and B components slowly 
ch anges, therefore an or bit of this pair was derived most recently 
by lAristidi et all (£999), who found a period of about 137 yr 
and a semi-major axis about 225 mas. The position angle of the 
pair has changed since this last paper by about 40°, hence a new 
analysis is required. Moreover, the combined analysis of the vi- 
sual orbit together with the times-of-minima variation due to the 
movement on this orbit and the radial-velocity variations can re- 
veal some of the other parameters of the orbit and also of both 
AB and C components. The more distant C and D components 
also belong to this multiple system, but these show no detectable 
mutual mo tion and can be a ssumed to be motionless. The MSC 
catalog by iTokovininl (I 1 9971) gives their periods 1 1 kyr and 591 
kyr, respectively. 

A distance to 65 U Ma was derived from the Hipparcos data. 
Perrv man et al.l (1 19971) gives the val ue of distance d - 2 46 + 
108 pc, while the new reduction of Ivan Leeuwe i (12001 pre- 
sented the value d = 212 + 30 pc. The c ombined spect r al typ e 
of the AB pair was classified as A3Vn bv lCowlev et alJ (II 969). 
while the spectrum o f the D component is A2p dSlettebald [l963). 
and lJoshi et al.l d2010l) carried out an analysis of this star. On the 
other hand, the spectral classification of the C component has 
never been performed. Therefore, the mass of the C component 
is also poorly constrained and the only information about this 
star that we have is a rough estimate of the magnitude difference 
(see below). 

3. Observations and data reduction 

In total, the target was observed on 82 nights: 29 nights for pho- 
tometry, and 53 nights for spectroscopy. The complete BVRI 
light curves of the eclipsing pair were obtained in 2010 at the 
private observatory of one of the authors (RU). However, owing 
to the relatively high brightness of the target, only a small 254- 
mm reflector of moderate focal length was used. This telescope 
was unable to separate the two 4" distant components, therefore 
the resultant LC was a composite AB+C light curve. The CCD 
photometric observations were obtained in standard B, V, and 
R filters according to the specification of lBesselH d 1990t) . All of 
the observations used for the LC were obtained with the same 
telescope and instrument setup, and the reduction was also iden- 
tical. Furthermore, the complete set of minima times used for the 
analysis is given in Appendix, two new minima were measured 
by Petr Svoboda, Czech Republic. 

The CCD spectra were obtained at Ondfejov observatory, 
Czech Republic, using the 2.0-m telescope equipped with a 
SITe-005 800 x 2000 CCD detector. These spectra cover a wave- 
length region 626 - 676 nm. All of them were secured be- 
tween March 2010 and July 2011 and have a resolving power 



Table 1. The parameters of the visual (A-B) orbit. 



Parameter 


Aristidi etal.(1999) 


This work 


Pa-b [yr] 


136.538 ± 8.4 


118.209 ± 0.690 


a A - B [mas] 


225 ± 18 


208.2 ± 9.7 


Aa-b [d] 




0.0428 ± 0.0023 


Ta-b 


2447140.9 ± 149.7 


2447516.9 + 126.8 


&a-b [deg] 


169.7 ± 4.6 


92.1 ±4.2 


<jJa-b [deg] 


26.9 ± 2.1 


202.7 ± 1.3 


i A -B [deg] 


39.7 ± 1.9 


38.1 ± 2.4 


eA-B 


0.531 ± 0.014 


0.504 ± 0.006 



R ~ 12700. Their S/N values range typically between 100 and 
300. 

For all of the spectra, the wavelength calibration was done 
using a ThAr comparison spectra obtained before and after the 
stellar spectra itself. The data reduction was performed fol- 
lowing the standard procedures of the data reduction package 
IRAFJ. The flatfields were taken in the beginning and end of 
each night and their means were used in the data reduction. 
After then, the radial velocities were obtained with the program 
SPEFO dHorn et all 1996llSkodal 19961) . using the zero point cor- 
rection by measuring the telluric lines. In total, 55 spectra were 
obtained in this way. 

All available data used for the analysis, the photometry, 
times of minima, spectroscopy, and the interferometric measure- 
ments are also listed (see the CDS tables). 

4. Visual orbit and the period analysis 

To begin with, we analyzed the visual orbit. The orbital motion 
influences the apparent period of the inner eclipsing pair, hence 
the periodic variation in the times of minima are analyzed ac- 
cording to the visual orbit parameters. 

Since its discovery as a double, 35 observations of the A- 
B pair have been obtained. These have been col lected in the 
Wash ington Double Star Catalog (hereafter WDSQ, Mason et~al1 
2001). We analyzed the data, obtaining the parameters of the vi- 
sual orbit presented in Table Q] and the final fit together with 
the data is given in Fig.Q] As one can see from Table [TJ the pa- 
rameters differ significantly in some aspects. Besides the higher 
precision, the most significant difference is found for the ori- 
entation of the orbit in space. It is obvious that the same fit to 
the data can be obtained with different sets of parameters when 
we only interchange the values of two parameters: (£2, lo) — » 
(Q + 180°, a) + 180°). However, when dealing with astrometric 
data set only, one cannot distinguish between these two identical 
solutions. The only way to do so is to use also the RV data, or 
the times-of-minima variation. 

For the minimum-times observations, we have only a limited 
set of data points. If we consider the period of the A-B pair to 
be about 118 yr, we have data for only about one-quarter of the 
orbital period covered with minima times at the present day. Yet, 
we can try to carry out an analysis of these data, by fixing the 
orbital parameters from astrometry (these in Table Q]). We have 
regularly observed the minima of this interesting target for the 
past four years to detect the period variation. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

2 http://ad.usno.navy.mil/wds/ 
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RA [arcsec] 

Fig. 1. Visual orbit of 65 UMa pair (A-B) as displayed on the sky. 
The eclipsing binary is placed in [0,0]. The dotted line represents 
the line of the apsides, while the dashed one is the line of the 
nodes. See section[4] 

Table 2. Final parameters of the short (Aa-Ab) orbit. 



Parameter 


Value 


PAa-Ab [d] 


641.5 ± 16.7 


A-Aa-Ah [d] 


0.00621 ± 0.00147 


T Aa-Ab 


2449615.4 ± 38.9 


^ Aa-Ab [deg] 


0.0 ± 15.2 


?Aa-Ab 


0.169 + 0.048 
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Fig. 2. Period variations of the eclipsing pair. Primary minima 
have been plotted as dots, and secondary as circles. The dashed 
line represents the 118 yr orbit, while the solid one is the 640 
day orbit. 



This led to an interesting finding that there is also an ad- 
ditional variation. We therefore analyzed our data set (43 data 
points in total) assuming two periodic terms. We used the Light- 
Time -Effect hypothesis (hereafter LITE, described e.g. bv llrwinl 
119591) . The results of our analysis can clearly be seen in Figs. [2] 
and[3] The long-term periodic modulation (blue) is caused by the 
118 yr visual orbit, while the short-period one is the newly dis- 
covered orbit, whose final parameters derived from our analysis 
are given in Table|2] This variation is clearly visible especially in 
the more precise recent data points after subtraction of the long- 
period term, see Fig. [3] Here we use the following labeling of the 
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Fig. 3. Period variation in the times of minima after subtraction 
of the 1 18 yr term. Only the variation caused by the component 
Ab and the most recent minima have been displayed. 



components: Aal and Aa2 for the eclipsing binary components, 
Ab for the 640 day orbit, and B for the 1 18 yr orbit (i.e. follow- 
ing the WDS notation). Thanks to the high precision of our new 
observations, the hypothesis of a non-circular orbit for 65 UMa 
eclipsing pair was ruled out. The complete list of times of min- 
ima together with the original BVRI photometry are given in 
electronic-only tables (see Appendix). There is a problem with 
some of the minima times, whose accuracy is not always given, 
hence we cannot perform a reliable chi-square test of the 640 
day hypothesis. However, using the weightening scheme for the 
data points, the LITE fit based on the 640 day hypothesis gives 
the sum of square residuals 0.00318, while disregarding the pos- 
sibility of a 640 day period the sum is 0.01222. 

Using the approach of combining the two LITE terms, one 
can also derive the parallax of the system independently of the 
Hipparcos value and the total mass of the system. The method 
is as follows: (Aa-b, oa-b) => n => M tot . To briefly describe the 
method, the amplitude of LITE and the angular semi-majo r axis 
of the visual orbit are directly connected via the parallax dMaveil 
Il990h . Using our new computed parallax and the Kepl er's third 
law, w e calculated the total mass of the system (e.g. iHilditch 
2001). The values presented in Table [2] and the LITE semi- 
amplitude Aa-b were calculated using this approach. This means 
that the values of pa-b, Ta-b, Qa-b, <^a-b, ia-b, and e A -B were 
fixed, but the parameters A A - B , PAa-Ab, A Aa -Ab, T Aa-Ab, ^ Aa-Ab, 
and eAa-Ab were fitted as free parameters. From this analysis, 
a new value of the parallax n = 4.28 + 0.49 mas was ob- 
tained, which yields the distance d - 234 + 29 pc. Such a value 
of parallax is slightly lower than the Hipparcos value (7THi P = 
4.72 + 0.58 mas). Using this new value of the parallax, we than 
computed the total mass of the system M tot = 8.25 + 1.85 M . 
lAristidi et al J (I 1 9991) found that M M =9.1 ± 1 1 .6 M . 

The relative motion of the component C around AB is very 
slow, but detectable. During more than 200 years of obser- 
vations, over 60 measurements were obtained (see the WDS) 
that revealed a change in position angle of about 5°. We ana- 
lyzed these data, determining a period of longer than 14000 yr. 
However, this result is very preliminary owing to the poor cov- 
erage of only 1/64 of the orbital period. 



5. Light and radial velocity curves 

To analyze the LC and RV curves, we had to consider a precise 
ephemeris of the eclipsing pair. These followed from the mini- 
mum times analysis and resulted in the elements for the primary 
minima 

HJD = 2455651.4491(5)+ 1.7304736(32) -E. (1) 
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Fig. 4. Radial velocity curves of the 65 UMa eclipsing pair. The 
upper plot shows the original RVs before correction for the 640 
day variation. The bottom plot shows the RVs of the eclipsing 
pair after the correction. 



xtuple system: 65 UMa 

Table 3. Parameters from the KOREL analysis. 



Parameter 


Value 


qAal-2 (=M Aa2 /M Aa i) 


0.995 ± 0.012 


K-a ci i [km/s] 


133.3 ±4.2 


K Aa2 [km/s] 


135.7 ± 4.2 


C/Aa-Ab (= M Ab 1 M Aa i +Aa2 ) 


0.69 ±0.11 


K Ab [km/s] 


19.9 ± 2.7 


qA-B (=M B /M A ) 


0.42 ±0.14 


K B [km/s] 


0.41 ± 0.30 
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Fig. 5. BVRI light curves of 65 UMa together with the final 
PHOEBE fit. The light curve parameters are given in Table [4] 



However, the analysis was not straightforward because of the ad- 
ditional components in the system. We assumed that the C com- 
ponent does not affect the spectra significantly, since its distance 
is about 4". The observing conditions and seeing were usually 
better than 2" during most nights. Hence, four other components 
were found to be present in ea ch of the spectra. We used the pro- 
gram KOREL dHadraval |2004|) to disentangle the spectra. 

To perform the spectral disentangling, the orbital elements 
of both orbits were fixed. Hence, the most crucial for the anal- 
ysis were the values of the mass ratios and amplitudes of the 
radial velocity curves. The ephemerides of the close eclipsing 
pair were also kept fixed because these are much more reliably 
known from the minima-times analysis. 

The parameters p Aa -Ab, T Aa - Ab , 0} Aa - Ab , e Aa - Ab , p A - B , T A - B , 
cj A -b, and e A s were fixed. The results of the other relevant pa- 
rameters using KOREL are listed in Table [3] Despite the results 
of the mass ratios not being very conclusive, we were able to 
make some preliminary estimations of the masses of the individ- 
ual components, as described below in Sect. [6] Since KOREL 
does not provide an error estimation, the errors in the individual 
parameters given in Table[3]resulted from the following analysis. 
Several solutions in KOREL were calculated, from which only 
those with^ 2 value closer than 5% from our best solution were 
considered. The errors in the parameters were assumed to be the 
maximum difference between these different solutions. 

The program KOREL enables us to obtain the RVs of the in- 
dividual components, which can be used for some further analy- 
sis. We used our knowledge of the ephemerides of the inner pair 
and the orbital parameters of the third body (i.e. 640 day orbit), 
to subtract the 640 day term from the RVs of the eclipsing pair. 
This can be clearly seen in Fig. [4] where the upper plot represents 
the original radial velocities, while the bottom plot represents the 
velocities after the subtraction of this term. We achieved signifi- 



cant improvement in the quality of the RVs, which could be used 
to perform a combined LC and RV analysis. 

The LC and RV curves of the e clipsing pair were an - 
alyzed using the program PHOEBE (Prsa&Zwitter, 2005), 
which is based on the Wilson-Devinney program (WD, 
IWilson & Devinnevlll971l) . The derived quantities are given in 
Table [4] while the LC has been plotted in Fig [5] The values of 
synchronicity F, for both eclipsing components were not derived 
from the combined LC and RV analysis. These were calculated 
from the spectra, and used to compute the values v sin i yielding 
the Fj values. The errors in F, are the standard deviations in F; 
measured for different spectra and different lines. The limb dark- 
ening was approximated using a linear cosine law, and the val- 
ues of Xi were interpolated from the tables given in Ivan Hammd 
d!993l) . In Table [4l we used the labeling of the two eclipsing 
components 1 and 2 in the indices instead of Aal and Aa2 for a 
better clarity. 

The primary temperature was fixed at valu e 7j_ = 800 K, 
which agrees with both the A7 spectral type dCoxl l2000h and 
the masses of the components (see Tabled). Therefore, it is le- 
gitimate to ask why a spectra l classification of an A3 star was 
made for the eclipsing binary dCowlev et al.U 19691) . This is due 
to presence of the component Ab on a 640 day orbit, whose spec- 
tral type i s probably of A3 and its l ight dominates the spectra. 
Moreover. iGarcia & Gimene9 d!986l) speculated that the eclips- 
ing binary components might be of A8-9 spectral type. The large 
value of the third light is also well-established owing to the three 
components that are present in the photometric aperture (in our 
notation, components Ab+B+C). 

The component Ab dominates the spectrum. This body has 
a well-defined orbit, hence its lines can also be plotted with the 
640 day period. On the other hand, Ab orbits around a bary center 
with the eclipsing pair. We can also plot the residuals from the 
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Table 4. Light curve parameters of 65 UMa eclipsing pair. 



Parameter 


Value 


Ti [K] 


8000* 


r 2 [K] 


7948 ± 20 


i [deg] 


86.5 ± 0.2 


Si = g2 


1.00* 


Ai =A 2 


1.00* 


Fx 


0.423 + 0.094 


Fi 


0.384 + 0.077 


U (B) [%] 


10.7 ± 0.4 


L 2 (B) [%] 


9.7 + 0.3 


L 3 (B) [%] 


79.6 ± 1.5 


U (V) [%] 


10.7 ± 0.5 


£2 (V) [%] 


9.8 + 0.5 


£3 (V) [%] 


79.5 ± 2.1 


£1 (R) [%] 


10.7 ± 0.7 


L 2 (R) [%] 


9.9 + 0.7 


L 3 (R) [%] 


79.5 ± 3.3 


U (I) [%] 


10.6 ± 0.4 


£2 (I) [%] 


9.9 + 0.4 


L 3 (I) [%] 


79.5 ± 1.8 


Derived quantities: 


Ry [Re] 


1.86 + 0.08 


R2 [Re] 


1.81 ±0.08 


Mi [M a ] 


1.74 ± 0.06 


M 2 [M Q ] 


1.71 + 0.06 


Note: * - fixed. 
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Fig. 6. Radial velocity curves on the 640 day orbit. The black 
squares stand for the Ab lines in the spectrum, while the red cir- 
cles represent the radial velocity residuals after subtraction of the 
eclipsing pair RV curve (filled for primary, open for secondary). 



LC fit from the eclipsing pair (from the original RVs), which 
should vary in anti-phase with respect to the Ab lines. This is 
shown in Fig|6] where we have plotted the 640 day fit to our 
spectra, which were acquired over two consecutive seasons 2010 
and 2011. 

Moreover, during the photometric monitoring of 65 UMa, 
two new variables were identified in the field . One of them 
was HD 103795 (spectrum K2III, according to Upgren 1962), 
w hile the other one was S AO 43913 (spectrum F0, according 
to ISlettebak & Stock1ll959l) . Neither was ever reported to be a 
variable, despite both having been observed by the Hipparcos 
satellite. However, our CCD photometry indicates that both are 
probably variable with amplitudes a few hundreds of magni- 
tude. SAO 43913 is probably a pulsating star (maybe 6 Set) 
with a period of about three hours, but the type of variability 
of HD 103795 remains unclear. 



sextuple system: 65 UMa 
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Fig. 7. Schematic structure of the whole system 65 UMa. 



6. Discussion and conclusions 

We have performed our first attempt to perform a detailed com- 
bined solution of all available data for 65 UMa, namely photom- 
etry, spectroscopy, and interferometry, obtaining quite a reliable 
picture of this unusual sextuple hierarchical system (see Fig. |7j. 

The inner close eclipsing pair consists of two almost- 
identical stars of A7 spectral type. This finding is consistent 
with the photometric indices B — V, V — R, and R - I being 
constant for the whole phase of the eclipsing binary at a level 
of 0.005 mag. The stars move on circular orbits with periods 
of about 1.73043 d, both being located on the main sequence. 
Thanks to the combined analysis, we were also able to compute 
its distance as d = 234 + 29 pc, independently of the Hipparcos 
satellite data. The 640 day orbit was confirmed by both the min- 
ima times and the RV variations. Applying the spectral disentan- 
gling and rough estimation of the mass ratios from this analy- 
sis, one can estimate the masses of the outer components. The 
Ab component is probably of Al spectral type and has a mass 
of about 2.4 + 0.4 M , from the total mass of the 118 yr vi- 
sual orbit, we can estimate the mass of the fourth component 
(in WDS named B), to be about 2.4 + 2.0 M Q . If we assume 
both these masses of Ab and B components, we can estimate 
the magnitude difference. I Aristidi et alj ( 1 1999b found this value 
to be 1.9 + 0.1 mag, while here we derive 0.7 ± 4.5 mag. The 
very large error is due to the large uncertainty in the mass of 
the fourth body. Another approach is to use the standard mass- 
luminosity relation and derive the individual luminosities of the 
components. Using this approach, we have plotted Fig. [8] where 
all components of the system 65 UMa are placed in the color- 
magnitude diagram. As one can see, the two eclipsing binaries 
are slightly under-luminous, while the D component seems to be 
over-luminous. The s ame finding about its higher luminosity wa s 
found elsewhere, e.g. lJoshi et alj d2010h or lAuriere et alj (2007). 

However, some properties of the Aa-Ab orbit remain un- 
clear, such as the inclination angle between the orbits. We can 
do some rough estimation of this quantity. The KOREL K^ a 
value and the predicted amplitude of radial-velocity variations 
from the LITE^^ are connected via sin iA a -Ab- Hence, we ob- 
tain i Aa-Ab ~ 47°, which lies well between the inclinations i of 
the eclipsing pair and the ia-b of the visual orbit. Nevertheless, 
its error is large but this is still only an estimation. We can also 
compute the predicted minimal angular separation of the Aa-Ab 
pair for a prospective interferometric detection. This resulted in 
about 1 1 mas, which is very favorable for modern stellar interfer- 
ometers, because the magnitude difference between the Aa and 
Ab components should also be rather low. On the other hand, the 
angular separation of the eclipsing pair components is still rather 
low, at about only 0. 18 mas. 
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Fig. 8. Color-magnitude diagram for all components of the sys- 
tem. Their position is compared with the Hipparcos stars (small 
black dots). 



Dealing with a multiple system, we should also consider the 
nodal period of the close pair and the 640 day orb it, hence the 
chang e in the inclination of the eclipsing binary (Soderhielm, 
fl975l) . The most crucial here is the ratio of periods pAa-Ab 2 IP, 
which implies that the nodal period was about 650 years, a du- 
ration that should be practical to observe. Unfortunately, we do 
not have a complete set of orbital parameters of the Aa-Ab pair, 
so this is only first rough estimation. However, this nodal period 
is not too long and potentially detectable. Further observations 
would help us to detect the change in the eclipse depths. Despite 
these being rather shallow, this effect was detected in only nine 
other cases, hence it would be interesting to reattempt detections, 
especially with the modern ultra-precise satellite photometry. 

Nevertheless, 65 UMa is a rather unusual system, we 
presently know of only 11 other sextuple systems (see 
Eggl eton & Tokovininl l2~008). The mass ratio of close to unity 
for the inner pair see ms to agree with some t heore tical models 
of star formation, e .g.lDe lgado-Donat e et al] d2003h . Moreover, 
some studies (e.g. Tokovinin & Smekhov 2002) indicate that 
about one-third of a ll multiples are higher-order systems. 
iGoodwin et al.1 (120071) discussed a finding that there is a dif- 
ference between the number of observed and expected higher- 
order multiples (quadruples and higher). Perhaps the discovery 
of other systems similar to 65 UMa would diminish this discrep- 
ancy. 
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Appendix A: Data sets used for the analysis A.2. Photometry 

A.1. Minima times 



Table A.l. Heliocentric minima times used for the analysis. 



HJD-2400000 


Error 


Tvpe 


Filter 


Reference 


43936.4851 




Prim 


B 


Gimenez & Ouesada (1982) 


43937.352 




Sec 


B 


Gimenez & Ouesada (1982) 


44249.683 




Prim 


B 


Gimenez & Ouesada (1982) 


44275.6355 




Prim 


B 


Gimenez & Ouesada (1982) 


44334.4700 




Prim 


B 


Gimenez & Ouesada (1982) 


44557.7034 




Prim 


B 


Gimenez & Ouesada (1982) 


44629.5218 




Sec 


B 


Gimenez & Ouesada (1982) 


48500.4500 




Sec 


Hp 


Perryman et al. (1997) 


49771.44535 




Prim 


B 


Drozdz (1997) 


51716.4305 


0.0004 


Prim 


BV 


Bulut & Demircan (2003) 


52381.7955 


0.0003 


Sec 


V 


Zasche et al. (2009) 


52381.7958 


0.0003 


Sec 


B 


Zasche et al. (2009) 


53399.2797 




Sec 


V 


Nagai (2006) 


53412.2596 




Prim 


V 


Nagai (2006) 


53499.6511 


0.0014 


Sec 


c 


Kraici (2006) 


53550.7037 


0.0017 


Prim 


c 


Krajci (2006) 


54166.7367 


0.0005 


Prim 


c 


Kraici (2007) 


54508.4870 


0.0012 


Sec 


R 


Zasche et al. (2009) 


54508.4886 


0.0005 


Sec 


BV 


Lampens et al. (2010) 


54514.5462 


0.0014 


Prim 


BV 


Lampens et al. (2010) 


54515.4132 


0.0013 


Sec 


V 


Lampens et al. (2010) 


54521.4630 


0.0009 


Prim 


R 


Zasche et al. (2009) 


54834.68161 


0.0012 


Prim 


R 


Zasche et al. (2009) 


54932.45785 


0.0006 


Sec 


R 


P.Svoboda - this paper 


54944.5637 


0.0014 


Sec 


V 


Lampens et al. (2010) 


55259.49539 


0.00054 


Sec 


BVR 


R.Uhlar - this paper 


55279.39706 


0.00048 


Prim 


R 


P.Svoboda - this paper 


55304.48805 


0.00050 


Sec 


BVRI 


R.Uhlar - this paper 


55535.51187 


0.00151 


Prim 


BVI 


R.Uhlar - this paper 


55561.46825 


0.00038 


Prim 


VRI 


R.Uhlar - this paper 


55613.37982 


0.0001 


Prim 


R 


Bratetal. (2011) 


55618.57558 


0.00092 


Prim 


C 


R.Uhlar - this paper 


55645.38457 


0.00146 


Sec 


c 


R.Uhlar - this paper 


55651.44796 


0.00030 


Prim 


I 


R.Uhlar - this paper 


55664.42595 


0.00022 


Sec 


I 


R.Uhlar - this paper 


55683.46083 


0.00089 


Sec 


1 


R.Uhlar - this paper 


55689.51398 


0.00194 


Prim 


I 


R.Uhlar - this paper 


55881.58599 


0.00116 


Prim 


c 


R.Uhlar - this paper 


55913.59785 


0.00049 


Sec 


c 


R.Uhlar - this paper 


55953.40014 


0.00037 


Sec 


R 


R.Uhlar - this paper 


55978.49541 


0.00061 


Prim 


V 


R.Uhlar - this paper 


55984.55331 


0.00019 


Sec 


I 


R.Uhlar - this paper 


56004.45281 


0.00089 


Prim 


c 


R.Uhlar - this paper 
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A3. Spectroscopy A.4. Interferometry 



Table A.2. Radial velocities used for the analysis. 



HJD-2400000 


R Vl,orig 


RV 2 , 


orig 


RV 3 


RVi, 


zorr 


RV Z , 


■orr 


55700.4062 


-143 


91 


119 


58 


11.32 


-133 


04 


130 


45 


55692.4050 


101 


22 


-133 


27 


20.87 


112 


05 


-122 


44 


55689.3063 


85 


21 


-108 


39 


12.67 


96 


02 


-97 


59 


55690.4499 


21 


06 


-45 


60 


16.34 


31 


88 


-34 


78 


55688.3121 


-140 


57 


117 


50 


12.38 


-129 


78 


128 


30 


55683.5066 


14 


74 


-36 


74 


13.83 


25 


48 


-26 


00 


55672.3907 


-94 


32 


70 


58 


13.40 


-83 


76 


81 


15 


55671.3706 


7 


48 


-25 


50 


12.43 


18 


02 


-14 


96 


55671.6085 


96 


10 


-121 


61 


12.50 


106 


65 


-111 


06 


55661.4470 


115 


79 


-143 


33 


15.66 


126 


10 


-133 


01 


55664.4727 


8 


33 


-35 


73 


18.82 


18 


72 


-25 


34 


55664.5012 


22 


05 


-46 


35 


15.87 


32 


44 


-35 


96 


55644.5210 


-2 


91 


-16 


85 


12.67 


6 


87 


-7 


06 


55642.4545 


113 


53 


-139 


80 


15.78 


123 


24 


-130 


09 


55623.4771 


95 


38 


-129 


19 


20.87 


104 


25 


-120 


31 


55622.2782 


-121 


65 


89 


66 


18.61 


-112 


84 


98 


47 


55622.5849 


-133 


41 


102 


54 


21.31 


-124 


58 


111 


37 


55619.5604 


55 


86 


-54 


86 


1.34 


64 


53 


-46 


19 


55619.6977 


98 


67 


-119 


32 


8.98 


107 


35 


-110 


65 


55618.5445 


17 


42 


-23 


71 


8.33 


26 


04 


-15 


10 


55616.5299 


107 


53 


-132 


28 


14.95 


116 


03 


-123 


78 


55602.5732 


121 


49 


-147 


33 


15.85 


129 


16 


-139 


66 


55425.3137 


-101 


14 


133 


75 


-20.57 


-115 


30 


119 


58 


55405.3360 


146 


75 


-115 


08 


-23.82 


130 


75 


-131 


08 


55399.3473 


-104 


42 


142 


16 


-28.06 


-120 


75 


125 


82 


55399.3633 


-99 


44 


141 


93 


-27.44 


-115 


77 


125 


59 


55398.3451 


150 


67 


-115 


91 


-26.18 


134 


29 


-132 


29 


55392.3562 


-114 


19 


152 


26 


-25.77 


-130 


79 


135 


67 


55391.3411 


138 


33 


-109 


01 


-23.67 


121 


71 


-125 


62 


55386.3829 


134 


08 


-106 


71 


-22.72 


117 


40 


-123 


39 


55385.3757 


-122 


49 


143 


68 


-13.58 


-139 


18 


126 


99 


55381.3681 


72 


64 


-47 


05 


-18.85 


55 


96 


-63 


73 


55372.4806 


143 


17 


-107 


56 


-23.15 


126 


71 


-124 


02 


55364.4965 


-98 


38 


135 


79 


-18.84 


-114 


43 


119 


74 


55358.3523 


108 


38 


-82 


47 


-16.83 


92 


76 


-98 


08 


55316.4097 


-86 


88 


95 


64 


-1.42 


-97 


36 


85 


17 


55316.4322 


-70 


94 


87 


93 


-14.14 


-81 


43 


77 


45 


55296.3370 


122 


71 


-127 


49 


3.33 


115 


34 


-134 


86 


55295.3948 


-124 


67 


147 


41 


-14.47 


-131 


89 


140 


18 


55293.3630 


-51 


44 


57 


58 


-3.28 


-58 


35 


50 


66 


55293.3797 


-55 


91 


72 


52 


-9.43 


-62 


82 


65 


61 


55293.3993 


-64 


32 


83 


21 


-12.32 


-71 


23 


76 


30 


55293.4198 


-72 


79 


89 


54 


-7.76 


-79 


71 


82 


62 


55290.4201 


-91 


57 


105 


50 


-8.15 


-98 


02 


99 


04 


55290.4826 


-61 


63 


87 


50 


-15.13 


-68 


09 


81 


04 


55290.4973 


-48 


36 


82 


18 


-17.47 


-54 


83 


75 


71 


55290.5149 


-51 


31 


69 


63 


-10.02 


-57 


78 


63 


16 


55290.5316 


-36 


63 


66 


50 


-16.93 


-43 


10 


60 


03 


55290.6121 


-9 


17 


26 


90 


-12.01 


-15 


65 


20 


41 


55289.4655 


131 


18 


-110 


46 


-10.99 


124 


88 


-116 


76 


55289.5066 


118 


55 


-101 


70 


-8.47 


112 


24 


-108 


01 


55280.4506 


95 


34 


-80 


39 


-9.83 


90 


42 


-85 


30 


55280.6528 


143 


80 


-120 


04 


-17.29 


138 


85 


-124 


99 


55274.3815 


-73 


08 


89 


51 


-7.98 


-77 


08 


85 


51 



Table A.3. Positional measurements of A-B pair used for the 



analysis. 


Year 




p [arcs] 


<f> [de; 


l\ 


Reference 


1908 


31 


0.26 


2 


9 


Aitken (19081 


1919 


63 


0.32 


11 


2 


Aitken & Doolittle 09321 


1934 


25 


0.30 


23 


4 


Kuiper (19611 


1934 


43 


0.32 


22 





Aitken & Moore (1937) 


1945 


36 


0.24 


34 


7 


van Biesbroeck (1954) 


1958 


09 


0.22 


54 


6 


van den Bos (1960) 


1958 


38 


0.21 


57 





van den Bos (1960) 


1958 


42 


0.24 


57 


6 


van den Bos (1960) 


1962 


14 


0.18 


65 


2 


van den Bos (1962) 


1963 


546 


0.19 


70 


3 


Worlev (1971) 


1964 


150 


0.21 


69 





Worley, (1971) 


1981 


4619 


0.121 


142 





McAlister et al. (1984) 


1983 


4141 


0.111 


157 


8 


McAlister et al. (1987) 


1984 


0530 


0.110 


162 


7 


McAlister et al. (1987) 


1984 


3751 


0.112 


164 


6 


McAlister et al. (1987) 


1984 


3834 


0.108 


164 


4 


McAlister et al. (1987) 


1985 


4894 


0.112 


168 


3 


McAlister et al. (1987) 


1986 


4066 


0.103 


179 


6 


McAlister et al. (1989) 


1987 


2695 


0.102 


186 


3 


McAlister et al. (1989) 


1988 


1653 


0.097 


193 


4 


McAlister et al. (1993) 


1988 


2523 


0.095 


196 


1 


McAlister et al. (1989) 


1989 


1561 


0.105 


198 


7 


Fuetal. (1997) 


1989 


2272 


0.099 


201 


9 


McAlister et al. (1990) 


1990 


2701 


0.095 


212 


2 


Hartkopf et al. (1992) 


1990 


2756 


0.097 


212 


4 


Hartkopf et al. (1992) 


1991 


3213 


0.101 


223 


4 


Hartkopf et al. (1994) 


1991 


3269 


0.100 


221 


7 


Hartkopf et al. Q994) 


1991 


9053 


0.088 


227 


8 


Hartkopf et al. (1994) 


1992 


3071 


0.091 


231 


6 


Hartkopf et al. (1994) 


1993 


1970 


0.093 


241 


9 


Hartkopf et al. (1994) 


1997 


074 


0.121 


276 





Aristidi et al. (1999) 


1998 


430 


0.120 


286 





Prieur et al. (2001) 


2007 


3114 


0.181 


328 


5 


Hartkopf & Mason (2009) 



